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We argue that the resonant peak observed in neutron 
scattering experiments on superconducting cuprates and the 
peak/dip/hump features observed in ARPES measurements 
are byproducts of the same physical phenomenon. We ar- 
gue that both are due to feedback effects on the damping 
of spin fluctuations in a d— wave superconductor. We con- 
sider the spin-fermion model at strong coupling, solve a set of 
coupled integral equations for fermionic and bosonic propa- 
gators and show that the dynamical spin susceptibility below 
Tc possesses the resonance peak at fires oc ^^^. The scatter- 
ing of these magnetic excitations by electrons gives rise to a 
peak/dip/hump behavior of the electronic spectral function, 
the peak-dip separation is exactly flrea- 
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One of the most intriguing recent developments in the 
physics of high Tc materials is the realization that not 
only the normal but also the superconducting state of 
cuprates is not described by a weak coupling theory. In 
particular, ARPES experiments on Bi2212 have demon- 
strated ||l|,|[ that even in slightly overdoped cuprates at 
T <^ Tc, the spectral function A(k, uj) near (0, tt) does 
not possess a single quasiparticle peak at lu = \/A^ + e^, 
where Ak is the superconducting gap and et is the 
fermionic dispersion. Instead, it displays a sharp peak 
which virtually does not disperse with fc, a dip at fre- 
quencies right above the peak, and then a broad maxi- 
mum (hump) which disperses with k and gradually recov- 
ers the normal state dispersion ||l|. Simultaneously, the 
neutron scattering experiments on near optimally doped 
YBCO 111 and Bi22U [g at T < T^ have detected a 
sharp resonance peak in the dynamical structure factor 
S{q,fl) oc x"{Qi^) centered at q — Q — (tt, tt) and at 
frequencies ~ 40 meV. 

In this communication we show that the resonance 
peak in S{Q,fl) and the peak/dip/hump features in 
^(^,0;) can be explained simultaneously by strong inter- 
action between electrons and their collective spin degrees 
of freedom which near the antiferromagnetic instability 
are peaked at or near Q. Specifically, we demonstrate 
that a d-wave superconductor possesses propagating col- 
lective spin excitations at frequencies smaller than twice 
the maximum value of the d— wave gap. The propagat- 
ing spin modes give rise to a sharp peak in S{Q,Q) at 
Q. = fires OC ^^^ where f is the spin correlation length. 
The interaction with collective spin excitations yields the 
fermionic self-energy H^ which at T = has no imaginary 



part up to a frequency loq which exceeds the measured 
superconducting gap by exactly Qres- 

The point of departure for our analysis is the spin- 
fermion model for cuprates which is argued [p[ to be the 
low-energy theory for Hubbard-type lattice fermion mod- 
els. The model is described by 



7^ = ^ VF(k - kf^)4^Ck,„ + J2 XO '(q)SqS-q 
k.Q q 

9 Yl 4+q,a '^",/3 Ck,/3 • S_q . 
q,k,CK,/3 
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Here cj^ ^ is the fermionic creation operator for an elec- 
tron with crystal momentum k and spin a, at are the 
Pauli matrices, and g is the coupling constant which mea- 
sures the strength of the interaction between fermions 
and the collective bosonic spin degrees of freedom. The 
latter are described by Sq and are characterized by a bare 
spin susceptibility xo(q) = Xo^V(l + il - Qf(^)- 

Eq. (nh gives rise to fermionic and bosonic self-energies 
and is particularly relevant for fermions near hot spots 
- the points at the Fermi surface separated by Q. In 
cuprates, the hot spots are located near (0, tt) and sym- 
metry related points. The presence of hot spots is essen- 
tial for our consideration because the fermions near these 
points are mostly affected by the interaction with anti- 
ferromagnetic spin fluctuations, and at the same time, 
they produce the dynamical part of the spin propagator 
because a spin fluctuation with a momentum near Q can 
only decay into fermions near hot spots. 

The normal state properties of the spin-fermion model 
have recently been analysed and compared with the ex- 
periments [|6|J3] . It was argued that the experimental sit- 
uation in cuprates corresponds to a strong coupling limit 
R = g/vpC^^ ^ 1, where g — g^Xo is the measurable 
effective coupling constant. The clearest experimental 
indication for this is the absence of the sharp quasiparti- 
cle peak in the normal state ARPES data for optimally 
doped and and underdoped cuprates |^,0] . At strong cou- 
pling, a conventional perturbation theory does not work, 
but it turns out that one can single out the most diver- 
gent diagrams, and incorporate them into a new "mean- 
field" ground state. This "mean-field" theory becomes 
exact if one formally sets the number of hot spots A^ = 8 
to infinity. At finite N, the expansion around the new 
vacuum holds in in {l/N)logR, but the prefactors are 
very small p] such that for practical purposes, one can 
restrict with the N = oo theory except very near the 
antiferromagnetic transition. 



This N — oo theory has some similarities with the 
mean- field d = oo theories |g] : it incorporates the domi- 
nant (~ R) self-energy correction which depends only on 
frequency, and also includes the dominant bosonic self- 
energy which is the spin polarization bubble made of the 
renormalized fermions. The corresponding set of self- 
consistent equations is presented in Eq. (g) for a super- 
conducting state. The normal state results are obtained 
by setting A = 0. 

The key physical effect which the TV = oo theory de- 
scribes is the progressive desctruction, with increasing 
R, of the coherent quasiparticle peak. At the same time, 
the fermionic incoherence has no feedback on spin sus- 
ceptibility which still has a simple relaxational form: 
X-Hq,n) = xo^V(l + iq- Q?e - *«M/) where 
iOsf = {2tt/N) {vf(,~^)/R @. The absence of the feed- 
back effect on spins is a quite general consequence of the 
fact that fermionic self-energy, abeit strong, has no de- 
pendence on the quasiparticle momentum ||9|] . 

In the superconducting state, this argument does not 
apply any more because superconducting and normal 
state Green's functions have different momentum depen- 
dences As a result, the feedback effect on spins is present, 
and already at the N = oo level one has to solve a set 
of coupled integral equations for the fermionic propa- 
gator and the spin polarization operator. This is the 
key intent of the present work. We however will not at- 
tempt to self-consistcntly find also the pairing suscepti- 
bility which in the spin-fermion model results from mul- 
tiple spin-fluctuation exchanges in the particle-particle 
channel pffl. Instead, we assume that below Tc the pair- 
ing susceptibility is a conventional (5-function of a to- 
tal momentum and frequency of a pair with the d— wave 
amplitude A^.. In other words, we will not distinguish 
between the true superconducting gap and the pseudo- 
gap. The full consideration should indeed include pair- 
ing fluctuations into the self-consistent procedure. We 
will also neglect the processes which scatter fermions 
near (0, tt) into fermions with momenta along zone di- 
agonal where the d— wave gap is absent. The contribu- 
tions from these processes soften sharp features associ- 
ated with the fc— independent gap, but are likely to be 
small numerically as they involve high energy spin fluc- 
tuations with momenta far from Q. Still, however, we 
will fully explore the fact that for d^2_y2 symmetry of 
the gap, Afe+Q = -Afc. 

A simple experimentation shows that 
for a (5— functional pairing susceptibility, the N — oo 
theory in the superconducting state has the same self- 
consistent structure as in the normal state but operates 
with normal and anomalous fermionic Green's functions 
G(fc,cj) = Goik,uj)/{l + Al Go{k,uj) Go{-k,-Lo)) and 
F{k,Lu) = -iAk/il + Al Go{k,u;) Go{-k,-uj)) where 
GQ^{k,uj) — w— E^j— Vi?(k— ki?) is the zero-order Green's 
function for 1/A^ expansion. This Go{k,uj) contains a 
self-energy which results from an exchange of a spin fluc- 



tuation with x(g, ^) = XoCV(l + (9-Q)^C^-no) where 
the spin polarization bubble 11^ is by itself a convolution 
of GG and FF. This construction yields a set of two 
coupled integral equations 



Scj = 3ig 



2 / d^qdfl 



G{k + q,uj + Q)xiq,^) 



nn = -2mg^^ I --^ {Gik,Lo)Gik + Q,u + n) 



(27r)3 
+F{k,uj) F{k + Q,uj)) 
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These equations has to be solved to leading order in 1/A^ 
(this, e.g., eliminates fc— dependence in S). Higer order 
self-energy and vertex correcions are small in 1/A^ and 
we neglect them. 

It is instructive to consider first the solution of Eqs 
(g) in the weak coupling limit A <S^ g <^ vpS,^^- To 
first approximation, IIq can then be evaluated with the 
free fermion Green's functions. This has been done be- 
fore |ll| and we just quote the result: in the supercon- 
ducting state, IIq has both real and imaginary parts. The 
/mlln = for fl < 2 A, it jumps a.t fl = 2 A to nA/uJsf, 
and then increases and approaches Q/tOsf at 17 3> 2A. 
This behavior is similar to that in an s— wave supercon- 
ductor except for the jump which is absent in s— wave 
case and is directly related to the fact that A^+q ~ — A^. 

The Relifi can be obtained either directly or using 
the Kramers-Kronig relation. At fJ <C A, we have 
Ren{ui) = (tt/S) f]2/(Aws/). It diverges at 2A as 
Ho = (A/uJsf) log(2A/|2A — ri|) because of the jump in 
/mllo, and decreases at larger frequencies. Due to the 
divergence, Rellci reaches 1 at a frequency fires which 
is less than 2A, i.e., when /mlln is still zero. Explic- 
itly, n-rcs = 2A(1-Z) where Z ex g-'^^/Zt^A), Near f^^es, 
xiQ,^) oc Z/(yi— fires — iS)i i-e., the dynamical structure 
factor has a resonance peak. 

Consider next the fermionic spectral function. With- 
out self-energy corrections, A{k,u!) near a Fermi surface 
resonates at LUres — A. The self-energy gives rise to a 
fermionic decay. For an s— wave superconductor, the on- 
set frequency for a decay is 3A, and ImTi(k,u)) emerges 
as {uj — 3A)^/^ 1^. The presence of the resonance mode 
qualitatively changes this picture because a fermion can 
decay into this mode starting from coq < 3A. A sim- 
ple power counting shows that this process yields a fi- 
nite jump of ImYi^ at the onset frequency and hence the 
logarithmical singularity in ReT,^. The latter in turn 
increases the self-energy at uj '^ A and shifts down- 
wards uJres (which is the measured gap), and the on- 
set frequency for ImllQ which, as one can easily demon- 
strate, exactly equals 2u!res- The amounts of the shifts 
and the amplitude of the jump in ImTi^ can be obtained 
explicitly from Eqs. (0). We found ljq = 3A(1 — e), 
uores = A(l — e), and 5[ImY^^^^) = (7rA/log2) e, where 
e = (31og2VA^/(64V^) {g/Af/^ ^-^.f/{2A)_ 



We see that the d— wave form of the gap yields qual- 
itative changes in the system behavior compared to the 
s—wave case, but at small coupling these changes are 
exponentially small and can hardly be measured. In par- 
ticular, the resonance peak in S{Q,uj) should be smeared 
out already by a small experimental resolution. The weak 
coupling results are shown in Fig M 
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FIG. 1. The T = weak coupling behavior of the dynam- 
ical structure factor and the spectral function. Solid lines 
are schematic solutions of Eqs. (Q) broadened by experimen- 
tal resolution. Without resolution, the peaks are (5-functions 
as indicated by arrows. The insets show the spin polariza- 
tion operator Iln and the fermionic self-energy E^, (solid lines 
- imaginary parts, dashed lines - real parts). The vertical 
dashed lines denote logarithmical singularities. 

We now turn to strong coupling limit ^ 3> A 3> 
(ujsf gY^^- Solving Eq. (0) as explained below, we found 
that the peak frequency in A{kF,i^) (i.e., the measured 
gap) is now w^es = A '^ A^/^. The nonzero ImTlQ and 
ImT,^ appear respectively at 2 A and loq — A(l-l-a) where 
a ex {iVsf/^y^^ oc ^~^. The amounts of jumps in Imlln 
and ImT,^ both scale as a^^^ and completely disappear at 
^ = oo when LOsf = 0. Above the threshold, ImT,^ first 
increases as Im'Ei^ ex {ui — luoY where i^ — \/3 — 1, and 
then recovers the normal state, ^/lu behavior. Substitut- 
ing this S into G{uj)^ we found that it possesses a peak at 
LJres , a dip at Wo and a hump at ujhump = i^res (1+&) where 
b ~ (5/A)3/(4'') ^1/(21.) logi? oc ^-1/(2^) log^. At C = oo, 
peak, dip and hump positions coincide with each other, 
and the peak/dip/hump structures transformes into the 
edge singularty: A{uj) oc (w — ujres)~'^- 

Further, the fact that ImllQ = up to 2uJres ini- 
plies, via Kramers-Kronig relation, that at small frequen- 
cies ReUn oc ^^^/(ws/A). Substituting this result into 
S{q,uj), we find that it possesses a resonance peak at 
fires ^ (w5/A)i/2 r^ ^-1 < 2A. At q ^ Q, the peak 
disperses with q a.s il'^ — (jj^g^((I -I- ((g — Q)^)^, just as 
a conventional spin wave, untill ft reaches 2aJres, and at 



larger frequencies disappears due to damping. 
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FIG. 2. Same as in Fig |l| but at strong coupling. The 
resonance and onset frequencies are presented in the text. 
The spin resonance frequency Q,res oc (,~^ , is equal to the 
distance between the measured gap A and the dip frequency 
Wo- The hump frequency differes from A roughly by ^ ' . 

The strong-coupling behavior of S{Q,fl) and A{k,uj) 
is presented in Fig. g. We see that (i) S{q,uj) pos- 
sesses a sharp resonance peak at fires ^ £,~^ which shifts 
with underdoping to lower frequencies, and (ii) yl(fci?,w) 
possesses a quasiparticle peak at uj — LOres^ a dip at 
ldq = LOres + fires, whcre ImTi^^ first appears, and a broad 
maxumum at a somewhat higher frequency ujhump- As 
the momentum moves away from the Fermi surface, the 
spectral function for frequencies larger than ujq disperses 
with k and recovers the normal state, non-Fermi liquid 
form with a broad maximum at w '--^ £fc/5- The quasi- 
particle peak however cannot move further than ujq be- 
cause of a strong fermionic damping above threshold. We 
found that it gets pinned at wq and just gradually looses 
its strength with increasing k — kp. 

We emphasize that although the resonance frequency 
in S{Q, uj) continuously evolves from weak to strong cou- 
pling, the physics changes qualitatively between the two 
limits. At weak coupling, the peak is solely due to a jump 
in ImliQ. At strong coupling, the jump is almost gone, 
and the existence of peak is due to fl^ behavior of Reliii 
which is related to vanishing /mllQ below 2ojres- 

We now briefly discuss how we obtained these results. 
We first integrated over momentum in (^ and for R = 
g/{vF£.~^) ^ 1 obtained 
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We then performed self-consistent calculations: we first 
obtained regular terms in i?eE and Rell by expand- 
ing (^ ^ in the external frequencies. We obtained 
i?eS,^ ex w(g/A)i/2 and ReUn ex OV(Aws/). In both 
cases, the integrals are confined to frequencies ~ A where 
the system interpolates between normal state and super- 
conducting behavior, and for estimates, one can use in 
the integrands the known normal state results for T,^ and 
n^^. We then assumed that at some finite frequency ojq 
ImTiaj jumps from to some finite value, considered the 
onset frequency and the amount of the jump as input pa- 
rameters, and used Kramers-Kronig relation to calculate 
the logarithmicaly singular term in ReT,^ . Adding it to a 
regular ReYt^j oc w, we find tOres where lo + ReT,^^^^^ = A. 
Substituting next ReY,^ into (Q) and using the spectral 
representation for Jmllo , we find the threshold frequency 
for /mlln at 2f2res and the amount of the jump at the 
threshold. We then again use Kramers-Kronig relation to 
calculate a logarithmically singular contribution to i^ell, 
add it to a regular Reliu oc il'^, and subsitute the result 
into (^ for Y,^. Using again the spectral representa- 
tion for ImT,^, we find two self-consistent equations for 
threshold frequency luq and for the amount of the jump at 
the threshold. These calculations are not exact because 
we only know regular ReY^j and -Rello upto overall nu- 
merical factors, but they yield correct functional forms 
of the resonance and threshold frequencies. 

We now compare our results with the data. Qualita- 
tively, the peak/dip/hump behavior of A{k,Lu), the ab- 
sence of the peak dispersion with k, and the presence 
of the dispersing resonance mode in S{q, lo) all agree 
with the ARPES and neutron measurements in YBCO 
and Bi2212 [|l|-0|. More quantitatively, we predict that 
the peak-dip separation in A{k, cu) at a hot spot exactly 
equals to the resonance frequency in S{Q,fl). Experi- 
mentally, in near optimally doped, Tc — 87K Bi2212, 
ujQ — ujres ~ A2meV |l|]. Recent neutron scattering 
data ID on Bi22l2 with nearly the same Tc = 91K 
yielded Qres — 'iSmeV, in full agreement with the theory. 
We also predict that the peak, dip and hump frequencies 
should tend to the same value as ^ ^- oo. At large but fi- 
nite ^, we predict that {oJhump/^res — 1) c>c (cjg/cj^gs — 1)^ 
where /3= (a/3 -h l)/4 ~ 0.7. 

We now connect our work to earlier studies. That the 
interaction with a nearly resonant collective mode peaked 
at Q explains the ARPES data has been known for some 
time, and qualitative arguments have been displayed first 
in IJI4] and then in m,. Rcf jl| also conjectured that the 
peak-dip separation may be related to a neutron peak 
frequency. It has been also realized earlier that in a 
d— wave superconducting Fermi gas, S{Q^ fi) contains 
a resonance peak exponentially close to 2 A |jll|. From 
this perspective, the key intension of this work was to 
present the quantitative results for cuprates by perform- 
ing actual strong coupling calculations, and to explicitly 
relate ARPES and neutron scattering data. 



Morr and Pines |^3[ obtained the spin-wave like disper- 
sion in x((7, f2) below Tc by phenomenologically adding 
the VP' term to the bare susceptibility. This term should 
be in the form Q? / ep as it can only come from fermions 
located far away from the Fermi surface. We have demon- 
strated that at T < Tc, the spin-fermion model of Eq ffl) 
by itself generates an fJ^/A term which for A ^ e^ com- 
pletely overshadows a possible bare fi^ term. 

Morr and one of us M considered an approximate solu- 
tion of Eq. (H) assuming that Ho still has the same purely 
relaxational form iil./u!sf as in the normal state, but Lu^f 
is momentum dependent. Comparing our results with 
||7|, we found that the approximate solution for A{k,Lu) 
captures the main features of the full solution, i.e., the 
peak/dip/hump structure of ^(^,0;), but yields incorrect 
values of the peak and dip frequencies for A ^ LUgf . 

J. Brinckman and P. Lee studied the evolution of the 
resonance peak within the slave boson theory [Q . Their 
philosophy and the results are similar to ours. 

To summarize, we considered the superconducting 
phase of cuprates and demonstrated that the reso- 
nance peak in the dynamic structure factor and the 
peak/dip/hump structure of the electronic spectral func- 
tion near (0, tt) can simultaneously be explained by a 
strong spin-fermion interaction. The peak-dip separa- 
tion at a hot spot exactly equals to the resonance neu- 
tron frequency and vanishes at ^ = cxd. The peak- hump 
separation also vanishes, but with a smaller power of ^~^. 

It is our pleasure to thank G. Blumberg, A. Finkelstein, 
R. Joynt, G. Kothar, A. Millis, D. Morr, M. Norman, D. 
Pines and J. Schmalian for useful conversations. The 
research was supported by NSF DMR-9629839. 



[1] 
[2] 
[3] 

[4] 
[5] 

[6] 
[7] 

[8] 

[9] 

[10] 

[11] 



M. R. Norman et al, Phys. Rev. Lett. 79, 3506 (1997). 

Z-X. Shen et al, Science 280, 259 (1998). 

H.F. Fong et al, Phys. Rev. B 54, 6708 (1996); P. Dai et 

al, Science 284, 1344 (1999). 

H.F. Fong et al, Nature 398, 588 (1999). 

P.B. Littlewood and CM. Varma, Phys. Rev. B 46, 405 

(1992). 

A. Chubukov, Europhys. Lett. 44, 655 (1997). 

A. Chubukov and D. Morr, Phys. Rev. Lett. 81, 4716 

(1998). 

A. Georges et al. Rev. Mod. Phys., 68, 13 (1996). 

L. Kadanoff, Phys. Rev. 132, 2073 (1963). 

P. Monthoux and D. Pines, Phys. Rev. B 47, 6069 (1993) 

D.J. Scalapino, Phys. Rep. 250, 329 (1995). 

D.Z. Liu, Y. Zha and K. Levin, Phys. Rev. Lett. 75, 

4130 (1995); I. Mazin and V. Yakovenko, ibid 75, 4134 

(1995); (1991); A. Millis and H. Monien, Phys. Rev. B 

54, 16172 (1996); N. Bulut and D. Scalapino, Phys. Rev. 

B 53, 5149 (1996). 



[12] Z-X. Shen and J.R. Schrieffer, Phys. Rev. Lett. 78, 1771 

(1997). 
[13] D.K. Morr and D. Pines, Phys. Rev. Lett. 81, 1086 

(1998). 
[14] J. Brinckmann and P.A. Lee, Phys. Rev. Lett. 82, 2915 

(1999). 



